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SPANNING TREES IN RANDOM SERIES-PARALLEL GRAPHS 


JULIA EHRENMULLER AND JUANJO RUE 


Abstract. By means of analytic techniques we show that the expected number of spanning 
trees in a connected labelled series-parallel graph on n vertices chosen uniformly at random 
satisfies an estimate of the form 

Se -"(l + o(l)), 

where s and g are computable constants, the values of which are approximately s ~ 0.09063 and 
g ps 2.08415. We obtain analogue results for subfamilies of series-parallel graphs including 
2-connected series-parallel graphs, 2-trees, and series-parallel graphs with fixed excess. 


1. Introduction 

The study of spanning trees and their enumeration is a central question in graph theory and in 
combinatorial optimization. It is well known that the number of spanning trees of a given graph G 
is an evaluation of its Tutte polynomial (see for instance El)- A lot of research has been devoted 
to study estimates of this number when dealing with restricted graph families. For instance, 
various results have been obtained for regular graphs and for graphs with degree constraints (see 

e.g. p EH EDI ED)- 

The enumeration of graphs with a distinguished spanning tree has also been studied extensively 
in the context of planar maps (namely, embedded connected planar graphs in the sphere, see 
Schaeffer’s Chapter at m for an introduction to this area). The first result in this context was 
obtained in the sixties by Mullin who was studying the number of rooted planar maps on n edges 
with a distinguished spanning tree [331 . Mullin determined that this number is C n C n + 1, where 
C n stands for the n-th Catalan number. Such formula was explained later by Cori, Dulucq and 
Viennot by means of Baxter permutations ESI, and by Bernardi using a direct bijection with 
pairs of plane trees with n and n + 1 edges, respectively (see 0)- Recently, Bousquet-Melou and 
Courtiel investigated the enumeration of regular planar maps carrying a distinguished spanning 
forest, as well as the connections of these counting formulas with statistical models as the Potts 
model m (see also [3] for the connection of spanning trees on maps and the Tutte polynomial). 

In this paper we study spanning trees in series-parallel graphs. A graph is series-parallel (or 
SP for short) if it is /q-minor free. Over the past few decades, SP graphs have been extensively 
studied from various points of view both in graph theory and computer science. In particular, 
being a subclass of planar graphs and a superclass of outerplanar graphs, SP graphs turned out 
to serve well as a pre-stage for the analysis and study of problems on planar graphs. Indeed, the 
family of SP graphs is the prototype of the so-called subcritical graph class family (see e.g. HUES]). 
In a typical connected graph in such a family, maximal 2-connected subgraphs (also called blocks) 
are small compared with the total size of the graph. This behaviour arises as a consequence of a 
subcritical composition phenomenon which appears in the specification of the generating functions 
associated with connected graphs of the family. 

In this paper we focus on enumerative problems defined on SP graphs. To this purpose, let us 
quickly state the following alternative definition of SP graphs that gives more insight into their 
structure and also justifies their name. Let G be a graph and let s and t be two of its vertices. We 
say G is series-parallel with terminals s and t if G can be turned into the single edge {s,t} by a 
sequence of the following operations: replacement of a pair of parallel edges (i.e. edges that have 
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Bordeaux in November 2014. J. R. was partially supported by the FP7-PEOPLE-2013-CIG project CountGraph 
(ref. 630749), the Spanish MICINN projects MTM2014-54745-P and MTM2014-56350-P, the DFG within the 
Research Training Group Methods for Discrete Structures (ref. GRK1408), and the Berlin Mathematical School. 
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two common endpoints) by a single edge, or replacement of a pair of series edges (i.e. non-parallel 
edges that have a common endpoint of degree 2) by a single edge. A graph G is 2-terminal series- 
parallel if there are vertices s and 1 in G such that G is series-parallel with terminals s and t. 
Finally, a graph G is series-parallel if and only if each of its 2-connected components is a 2-terminal 
series-parallel graph (see e.g. fl2|). 

Also, SP graphs are known to be the class of graphs of treewidth at most 2 (see e.g. El)- 
Edge-maximal SP graphs (i.e. graphs which cease to be SP whenever any missing edge is added) 
are exactly the class of 2-trees, which can be defined in the following way. A single edge is a 
2-tree. If T is not a single edge, then T is a 2-tree if and only if there exists a vertex v of degree 2 
such that its neighbours are adjacent and T — v is also a 2-tree. Conversely, every subgraph of 
a 2-tree is a SP graph. In particular, SP graphs are at most 2-connected since 2-connected SP 
graphs always contain a vertex of degree 2. 

From now on, unless stated otherwise, all graphs under study are labelled and simple. By a 
random object of a given family we mean an object chosen uniformly at random from all the 
elements of the same size, e.g. graphs on the same number of vertices. In the present paper we 
study enumerative properties of spanning trees and spanning forests on random SP graphs. Before 
stating our results, let us survey some relevant related investigations. 

One can easily verify that the number of edges of an n-vertex 2-tree is precisely 2 n — 3. In the 
same context, Moon [[32] showed that the number of 2-trees on n vertices is equal to (™) (2n —3) n_4 . 
The enumeration of SP graphs is, however, more involved. Bodirsky, Gimenez, Kang, and Noy 
proved in [8] that the number of connected SP graphs on n vertices is asymptotically of the form 

c s n~ 5/2 g~ n n \, 

where c s « 0.00679 and g s ss 0.11021 are computable constants. In the same paper they showed 
that the number of edges in a random connected SP graph is asymptotically normally distributed 
with mean asymptotically equal to nn and variance asymptotically equal to An, where n ~ 1.61673 
and A ~ 0.2112 are again computable constants. 

Building on these results, a lot of research has been done to understand the qualitative picture 
that emerges when studying a random SP graph with a fixed number of vertices. The maximum 
degree and the degree sequence of a random SP graph have been studied in [2D] and mm, 
respectively. Drmota and Noy m investigated several extremal parameters in subcritical graph 
classes, which include the class of SP graphs. They showed, for instance, that the expected 
diameter D n of a random connected SP graph on n vertices satisfies c\y/n < E[D„] < C 2 \/n log n 
for some positive integers c\ and C 2 . The precise asymptotic estimate has been proved very 
recently by Panagiotou, Stufler, and Weller [36] to be of order Q{^/n). In the same work, the 
authors exploited this fact to prove that in subcritical graph classes, and in particular in SP 
graphs, the normalized metric space {V (G ), do / y/n) (where dci;u,v) is the number of edges in a 
shortest path that contains u and v in G) converges with respect to the Gromov-Hausdorff metric 
to the Brownian Continuum Random Tree multiplied by a constant scaling factor that depends 
on the class under study (see [36]). 

Our results. In the present paper we study the number of spanning trees in a random (connected 
or 2-connected) SP graph on n vertices. In particular, our main result is a precise estimate of the 
expected number of spanning trees. 

Theorem 1.1. Let X n and Z n denote the number of spanning trees in a connected, respectively 
2-connected, labelled SP graph on n vertices chosen uniformly at random. Then, 

E[X n ] = sg~ n { 1 + o(l)), where s « 0.09063, g~ x w 2.08415, 

E [Z n }= pw~ n {l + o{l)), where p « 0.25975, w~ x « 2.25829. 

The previous analysis is done over all (connected, 2-connected) SP on a given number of vertices. 
However, we also address the study of extremal situations. First, we can also particularize the 
computation of the expectation in the case of a random 2-tree on n vertices, which maximizes 
the number of edges in an n-vertex SP graph. In this case, the expected value of the number of 
spanning trees is slightly bigger than the one in Theorem 0 
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Theorem 1.2. Let U n denote the number of spanning trees in a labelled 2-tree on n vertices 
chosen uniformly at random. Then, the expected value of U n is asymptotically equal to S 2 Qf n , 
where S 2 ~ 0.14307 and gf 1 ~ 2.55561. 

Finally, we study SP graphs with few edges. More precisely, we elaborate the expected number 
of spanning trees in a random connected SP graph on n vertices and fixed excess equal to k, where 
k is a integer that does not depend on n. Recall that the excess of a graph G is defined as the 
number of its edges minus the number of its vertices (by fixed we mean that it does not grow as 
a function of n). Our result is a polynomial estimate (in n ) of the expected number of spanning 
trees: 


Theorem 1.3. Let k > 1 be a fixed integer. Let X n j~ denote the number of spanning trees in a 
connected labelled SP graph, on n vertices and with fixed excess equal to k, chosen uniformly at 
random. 

Then, when n is large enough, 


E[A' n>fc ] 


m 


P(3fc/2) 

r(2fc + 1/2) 





where the function c(k) satisfies that, for k large 


c(k) = £7 fc (l + o(l)), 


( 1 ) 


with c ss 0.90959 and 7 1 « 2.60560. 


The previous formulas must be understood in the following way: we fix k and we let n tend 
to infinity. Additionally, if k is sufficiently large, we can get the approximation of c(k) stated in 
the second part of Theorem 1.3 Indeed, the term o(l) in Equation |l]) is polynomially small in k 

In order to deduce these expressions in Theorem |1.3| we analyse weighted cubic SP multigraphs 
on 2k vertices. These objects are reminiscent of the work |28] and are building on previous 
enumerative results on simple cubic planar graphs [5], see Section [ 6 ] for definitions and details. 
The asymptotic estimate stated in Theorem 1.3 arises when getting asymptotic estimates (in terms 
of k) for the number of such multigraphs. 


Organisation. The rest of the paper is organised as follows. In Section [2] we introduce the 
essential combinatorial and in Section [3] the essential analytic definitions, techniques, and results 
that we use. The proof of Theorem |l.l| is then presented in Subsection |4.1| of Section[4j In the same 
section we also analyse the behaviour of the growth constant of the expected number of spanning 
trees in a random connected SP graph if we fix its edge density (Subsection |4.2| ) and comment on 
the variance of the number of spanning trees in a random 2-connected SP graph (Subsection 4.3). 
Next, Section[5]is devoted to the analysis of 2-trees and the proof of Theore m|1.2| Then, Section 6 
deals with SP graphs with fixed excess and presents the proof of Theorem |1.3[ Finally, Section 7 
contains some concluding remarks and open problems. 


2. Combinatorial Preliminaries 

Notation. Our combinatorial and analytic notation is standard and follows [53] ■ In particular, 
given a formal power series of exponential type A(x) = (C u >o °' n X \ (EGF for short), we use the 
notation [ x n ]A(x ) to indicate the n-th coefficient of A{x). Given a bivariate function A{x,y), 
we denote the partial derivative of A(x,y) with respect to x and y by A x (x,y) and A y (x,y), 
respectively. However, we will usually use the notation A'(x,y) to denote A x (x,y). We write 
o n ~ b n whenever linin^oo a n /b n = 1. Throughout the paper log denotes the natural logarithm. 
In our setting, we use the variable x to mark vertices and the variable y to mark edges. These 
variables are exponential and ordinary, respectively. 
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Graph decompositions. The main ingredients in our proofs, from an enumerative combinatorial 
point of view, are the Symbolic Method, generating functions, connectivity decompositions, and 
an extension of the Dissymmetry Theorem to tree-decomposable classes. In this section we review 
the essential definitions and results related to these topics. For further details, in particular for an 
introduction to the Symbolic Method, we refer to the book Analytic Combinatorics by Flajolet 
and Sedge wick [23] - 

Let C be a class of connected graphs with the property that a graph is in C if and only if 
all its 2-connected and 3-connected components are in C. Observe that for instance the class of 
connected SP graphs carrying a distinguished spanning tree shares this property. Let c n ^ m denote 
the number of graphs in C with n vertices and m edges. The associated (mixed) exponential 
generating function (or EGF for short) is the formal power series 

C{x,y)= V C nm r,y m , 

' n\ 

m,n>0 

where x and y mark vertices and edges, respectively. 

Similarly, let denote the number of 2-connected graphs in C with n vertices and m edges 
and let B{x , y) be its associated EGF. A connected graph rooted at a vertex can be obtained from 
a set of rooted 2-connected graphs, where the root is not labelled and where every other vertex is 
substituted by a connected graph rooted at a vertex. Using the Symbolic Method, this translates 
into the following relation between C(x,y) and B{x,y) (see also [23]). 

xC'(x, y)=x exp (B'(xC'(x, y), y)) . (2) 

Following [40] . a network is defined as a simple graph with two distinguished vertices, that 
are called 0-pole and oo-pole and do not bear a label, such that adding an edge between the two 
poles creates a 2-connected multigraph. If there is an edge joining the two poles, it is called root 
edge. The root edge defines the two poles, which are usually denoted by 0 and oo (initial and final 
vertex of the root edge, respectively). Let D(x,y) denote the EGF associated with networks. The 
following equation, shown by Walsh in m, reflects the relation between B(x,y) and D{x,y). 

2(1 +y)B v (x,y) = x 2 {D{x,y) + 1). (3) 

The left-hand side in Equation ([3]) corresponds to the family 2-connected graph rooted at a directed 
edge that might not be present in the graph and the right-hand side corresponds to the family of 
networks (possibly empty) where in addition a label is given to the two poles. 

A trivial network consists of the two poles and of the root edge. Following the ideas of [38] , 
we further distinguish between three types of networks, namely series , parallel , and h-networks 
as follows. A series network S can be obtained from a directed cycle with a distinguished edge 
(which defines the two poles of the network) by replacing every other edge by a network, and 
finally removing the distinguished edge. A parallel network P arises from merging at least two 
non-trivial networks, the root edge of each of them being not present, at their common poles. 
In this family, the root edge joining the two poles of P might not be present in P. Finally, an 
//.-network is obtained from a 3-connected graph H rooted at an oriented edge, by replacing every 
edge of H apart from the root edge by a network. As in the parallel case, here the root edge might 
not be in the network. 

Trakhtenbrot [38] showed that a network is either trivial, series, parallel, or an //-network, and 
Walsh m translated this decomposition into counting formulas. In SP graphs, the set of h- 
networks is empty, so in the rest of the paper we deal only deal with series and parallel networks. 

Let us finally mention that our definition of network slightly differs from Trakhtenbrot’s. Indeed, 
in [38] series networks could contain the root edge. In our work, series networks containing the 
root edge (in Trakhtenbrot’s sense) are always considered to be parallel. This convention would 
arise to be helpful when dealing with spanning trees. 

In Section [4] we aim for a precise asymptotic estimate for the number of spanning trees in 
random SP graphs. For this purpose, we will enumerate the class of connected SP graphs with a 
distinguished spanning tree. The main idea is to give a complete analytic analysis of the gener¬ 
ating function associated with this class using the relations to the class of 2-connected SP graphs 
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and to the class of networks, both carrying a distinguished spanning tree. Using Equation § 
would imply integration steps that are known to get difficult when considering enriched classes 
of graphs. Chapuy, Fusy, Kang, and Shoilekova m found, however, a convenient combinatorial 
trick to forgo this integration step by using the dissymmetry theorem for tree-decomposable classes 


(Theorem 2.1l and by using the grammar for decomposing graphs into 3-connected components 
that they developed in m- Theorem [2d] will also serve us well in Section [5} 

A class A of graphs is tree-decomposable if for each graph G £ A we can define a tree t(G) 
associated with G. Let A 0 denote the class of graphs G in A where r(G) has a distinguished vertex. 
Similarly, denote by A 0 - 0 the class of all graphs G in A where t(G) carries a distinguished edge. 
Finally, let A a ->o be the class of all graphs G in A where an edge of r(G) is directed. The 
Dissymmetry Theorem for trees by Bergeron [2] allows to express the class of unrooted trees in 
terms of classes of trees with a distinguished vertex, edge or with a directed edge. This theorem 
can be extended to tree-decomposable classes in the following way (see e.g. ED- 

Theorem 2.1 (Dissymmetry Theorem for tree-decomposable classes). Let A be a tree-decompo¬ 
sable class of graphs. Then, 

A + Ao— K5 — Ao + Ao-O- 

Finally, let us briefly summarize Tutte’s decomposition |39) for decomposing 2-connected graphs 
into 3-connected components. For a thorough exposition we refer to |T4], 

Tutte’s decomposition is based on split operations and the structure obtained from this process 
is shown to be independent of the order of the operations. Roughly speaking, in every split 
operation we split the edge set of a graph G into two edge sets E\ and E 2 that only coincide in 
exactly two vertices, say u and v, and where G[E\\ is 2-connected and G[E 2 ] is connected modulo 
{u, v} (meaning that there exists no partition of E 2 into two nonempty sets E 2 and E'f such that 
G[E’ 2 ] and G[E'f] only intersect in u and v). Next we add a so-called virtual edge e between these 
two vertices. Then we split the graph along this virtual edge which yields two graphs G\ and G 2 
that correspond respectively to G[E\\ and G[E 2 ] with e now being a real edge. We say that G\ 
and G 2 are matched by the virtual edge e. 

The resulted structure is a collection of graphs that we call bricks. Tutte showed that there 
are only three types of bricks, namely ring graphs (R-bricks), multi-edge graphs (M-bricks), and 
3-connected graphs with at least 4 vertices (T-bricks). The class of ring graphs is defined as the 
class of cyclic chains of at least 3 edges and the class of multi-edge graphs as the class of graphs 
with exactly two labelled vertices that are connected by at least 3 edges. 

The RMT-tree of a graph G is defined as the graph t(G) the vertices of which are the bricks 
that result from Tutte’s decomposition applied to G. Two vertices in r(G) are connected, when 
the corresponding bricks are matched by a virtual edge. It was shown by Tutte [35] that r(G) is 
indeed a tree and there are no two adjacent R-bricks nor two adjacent M-bricks. 

Let B be the class of all 2-connected graphs with at least 3 vertices. We denote by Br, Bm , and 
Bt the classes of graphs G in B such that the RMT-tree associated with G carries a distinguished 
R-vertex, M-vertex, and T-vertex, respectively. Moreover, let Br-m denote the class of graphs G 
in B such that an edge between an R-vertex and an M-vertex in the RMT-tree associated with G 
is distinguished. The classes Br-t , Bm-t, and Bt-t are defined analogously. Finally, let Bt^t 
be the class of graphs G in B such that an edge between two T-vertices is directed. 


Using Theorem 2.1 B satisfies the following equation as shown in M- 

B ~ B r + Bm + Bt — Br-m — Br-t — Bm-t — Bt^t + Bt-t- (4) 

In our work we consider only SP graphs. In particular, a SP graph does not have 3-connected 
components. This implies that SP graphs do not contain T-bricks, and in hence RMT-trees do 
not have T-vertices. So Equation Q is simplified to 

B ~ Br + Bm — Br-m- (5) 


3. Analytic Background 


The proofs in this paper are based on singularity analysis of generating functions. In this 
section we introduce the necessary analytic background. For the sake of completeness, we state the 
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results that we use, in particular a simplified version of the Transfer Theorems (Theorem |3.1[ ) and 
a simplified version for the singularity analysis of systems of functional equations (Theorem |3.2[ ). 
For more details, we refer to the books Analytic Combinatorics by Flajolet and Sedgewick |2.'li 
and Random Trees by Drrnota El- 

Given a univariate exponential generating function 

A(z) = ^a„ — 

n> 0 

we would like to determine an asymptotic estimate of the sequence (a n ) n > 0 - Pringsheim’s Theorem 
(see e.g. [23]) assures that generating functions with radius of convergence g and non-negative 
Taylor coefficients have a singularity at g , in particular a positive real dominant singularity. As 
shown in [23 (Theorem IV.7), the exponential growth of the sequence (a n ) n >o is therefore dictated 
by the smallest positive singularity g of A(x) in the sense that 


\x n ]A{x) ~ Q(n)g ", 

where 0(n) grows subexponentially, i.e. limsup n _ >00 |0(n)| 1 /" = 1. The subexponential term 
0(n) results from the nature of this singularity. The so-called Transfer Theorems , developed by 
Flajolet and Odlyzko [22] . provide us a convenient way to determine the subexponential term of 
[x n ]A(x). In particular, Theorem 3.1 is a special case of the Transfer Theorems in [23]. For this, 
we need the definition of dented domains. Given R, £ > 0 with R > f and 0 < <f> < 7r/2, the 
domain dented at f (which we write as A^(<(>, R )) is defined as 

A c (0, R) = {z G C : |z| < R, z ± (, | Arg(z - C)| > 01- 


Theorem 3.1. Let a G R\Z and let A{x) be analytic in a domain A e (</>, R) dented at the 
smallest positive singularity g of A{x). If, as x —» g in A e {(j),R), 


A(x)^c(^l- J ^j 

then 

[x n ]A(x) = -^-n Q_ 1 p""(l + o(l)), 

T(a) 

where T(a;) is the classical Euler Gamma function defined as T(a;) = J 0 °° t x ~ 1 e~ t dt. 


In this paper, the singular expansion of a generating function A{x) in a domain dented at a 
singularity g is always of the form 


A(x) = A 0 + A\X + A 2 X 2 + ... + A 2k+ iX 2k+1 + 0(X 2k+2 ), 


where X = \Jl — x/g. The even powers of X, being analytic functions, do not contribute to the 
asymptotic of [x n ]f(x). 

If A\ = A 3 = ■ ■ • = A 2k _i = 0 and A 2k+ i 0, then the number (2fc + l)/2 is called the singular 


exponent. Then, by Theorem 3.1 we get that 

[x n ]A(x) r 


T(a) 


with c = A 2 fc+i and a = — (2fc + l)/2. When dj / 0 we say that A(x) has a square-root expansion. 

Let us now turn to the asymptotic analysis of systems of functional equations. The main refer¬ 
ence for this topic is the paper [TBj. We include here a shortened version (see Section 2.2.5. in im 
for the more general statement). Assume that yi{x), ... ,y k {x) are generating functions satisfying 
a functional system of equations. We define y = (yi(x ),..., y k (x)), and the system satisfied by y 
is denoted by y = F(x; y), where F = (Fj ,... ,F k ) is a vector of functions. The dependency graph 
G = ( V. , E) associated with the system y = F(t; y) is an oriented graph the vertex set of which is 
V = {yi,... ,y k } and yfy f j is in E if and only if 7 ^ 0. The latter condition indicates that there 
is a real dependence between Fj and yj. A dependency graph is said to be strongly connected if 
every pair of vertices can be linked by a directed path. Using this terminology, we can finally state 
the following result: 
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Theorem 3.2 (Systems of functional equations [T7] . simplified version). Consider the functional 
system of equations y = F(x; y) satisfying that each yi is analytic at x = 0. Additionally, we 
require that each component of F is an entire function with positive Taylor coefficients, that it 
is not linear in the components yi and depends on x. Finally, we assume that F( 0; y) = 0 and 
F(x; 0) ^ 0. Assume also that the associated dependency graph is strongly connected. Denote by 
Ik the k x k identity matrix and by Jac(T') the Jacobian matrix associated with F and with respect 
to variables j/i,... ,j/fc. Assume that the system of equations 

y=F(x;y), 0 = det (I k - Jac(-F)) 

has a unique solution (xq, y 0 ) in the region of analyticity of each component of F. Then there is a 
unique solution y of the initial system of equations such that the components of y have non-negative 
Taylor coefficients and a square-root expansion in a domain dented at x = Xq- 


In order to obtain asymptotic estimates we need to assure that the dominant singularity is 
unique in a dented domain. This condition is usually satisfied whenever the counting formula 
A(x) under consideration cannot be written in the form A(x) = x k f(x r ) for non-negative values 
k > 0 and r > 2. More precisely, we say that a generating function A[x) is aperiodic if there 
exists a non-negative integer no such that \x n }A{x) > 0 for n > no- Observe that checking the 
aperiodicity condition is straightforward whenever we know that for each number of vertices there 
exist graphs in the family under study. The generating functions we consider in the forthcoming 
section (which are defined by an implicit equation, or by means of Theorem 3.2) will satisfy 
the aperiodicity condition by obvious combinatorial reasons. This will imply uniqueness of the 
dominant singularity. See m for details. 


4. Spanning trees in series-parallel graphs 

In this section we present the proofs of Theorem 1 1.1 1 determining a precise asymptotic estimate 
for the expected number of spanning trees in random connected SP graphs. Additionally, we 
elaborate the expected number of spanning trees in a random 2-connected SP graph of a given 
edge density. Recall that all graphs in this paper are considered to be labelled. In order to 
count spanning trees in SP graphs, we are concerned with the enumeration of SP graphs carrying 
a distinguished spanning tree. For this, let c„ jm and now denote the number of connected, 
respectively 2-connected, SP graphs with a distinguished spanning tree and let C(x, y) and B(x , y) 
be their associated counting formula, where x and y mark again vertices and edges, respectively. 
Furthermore, let D denote the class of series-parallel networks carrying a distinguished spanning 
tree and let D(x,y) be its associated generating function. 

4.1. Expected number of spanning trees in random SP graphs. The first step in our 
proof of Theorem [hi] is the enumeration of SP networks that carry a distinguished spanning tree. 
For the purpose of counting spanning trees in SP networks, we need to introduce the following 
auxiliary class. Let V denote the class of SP networks that carry a distinguished spanning forest 
with two components, each of which contains one of the poles. Let D(x, y) be its associated EGF. 

Recall that a network is either trivial, series or parallel. By convention, we assume that networks 
with a root edge are parallel. Therefore, we define the following classes of networks. Let S 
and S denote the class of series networks that carry a distinguished spanning tree, respectively 
a distinguished spanning forest with two components each of which contains one of the poles. 
We denote their associated EGFs by S(x,y) and S(x,y), respectively. Similarly, let V and V 
denote the class of parallel networks and which carry a distinguished spanning tree, respectively 
a distinguished spanning forest with two components, each of which contains one of the poles. 
Observe that in both families the root edge might be present. Their associated EGFs are denoted 
by P{x,y) and P(x,y), respectively. For the sake of readability we may omit the parameters 
whenever they are clear from the context. 

We start with elaborating relations between D , D, S, S, P, and P in order to obtain a suitable 
system of equations. One can easily verify that 

D(x, y) = y + S(x,y) + P(x, y), 


( 6 ) 
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and 

D(x,y) = y + S(x,y) + P{x,y). (7) 

Note that in Equation Q the variable y on the right-hand side corresponds to a trivial network 
with a distinguished spanning tree, whereas in Equation |7| it corresponds to a trivial network 
with a distinguished spanning forest that consists of two components of size 1. 

Let us now analyse series networks. Observe that a series network N can be decomposed into 
at least two networks, where the 0-pole of the i-th network is identified with the oo-pole of the 
(i + l)-th network. Equivalently, N can be decomposed into an ordered sequence formed by a 
network N 1 that is not series and an arbitrary network N" that are joined by a series operation. 
If N £ S, then each of these two networks contains a distinguished spanning tree. Therefore, we 
have 

S{x, y) = (D(x, y) - S{x, y))xD{x, y) = (y + P(x, y))xD(x, y ). (8) 


V\S V 



Figure 1. Decomposition of N £ S. 

If N £ S, then either N' £ V\S and N" £ V, or N' £ D\S and N" £ V. This translates into 
the following equation: 

S{X, y) = (D(x, y) - S{x, y))xD(x , y) + (D(x, y) - S{x, y))xD(x , y) 

= (y + P(x,y))xD(x,y) + (y + P{x,y))xD(x,y). (9) 


V\S V T>\S V 



Figure 2. Decomposition of N £ S. 


Finally, we analyse parallel networks. A parallel network can be described as a set of at least 
one series network if the root edge is present, or of at least two series networks, otherwise. If 
N £ V, we need to distinguish between the case that the root edge is present and the case that 
it is not. In the second case, all series networks are in S except for one which is in S. If in the 
first case the root edge is in the distinguished spanning tree of N, then all series networks are in 
S. If, on the other hand, the root edge is not in the spanning tree, then exactly one of the series 
networks is in S and all other networks are in S. Thus, we get 

P(x,y) = y(exp(S(x,y)) - l) + y(S{x,y)exp(S(x,y))) + S(exp(S{x,y)) - l). (10) 


If N £ V, then N can be decomposed into the root edge if present and into series networks in S 
that are joined by a parallel operation. If the root edge is present, then there is at least one other 
network. In the other case, there must be at least two. This gives rise to the following equation: 


P{x,y) = (exp {S{x,y)) - S{x,y) - l) + y(exp(S{x,y)) - l). (11) 

Using formal manipulations, we get that the system of Equations defines the following 

implicit expression for D(x,y): 

xD 2 \ ( „ {y(l + xD)-(l + y)D)(2 + xD) 


D= y + (l + y) 


1 + xD 


exp — xD 


(2/(1 + xD) + (1 + y)xD 2 )(l + xD) 2 


( 12 ) 


In order to study the singular behaviour of all the previous generating functions we could apply 
Dromta-Lalley-Woods methodology for systems of functional equations (see e.g. [23J). However, 
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0 0 0 





Figure 3. Decomposition of N G V. 


o 



o 



Figure 4. Decomposition of N € V. 


as in this particular case we have an expression for D{x,y ) not depending on any other variables 
but x and y, we will analyse Equation (12) in order to get the singular behaviour of D[x , y). The 
following theorem is reminiscent to Lemma 3.3. in [25]: 


Lemma 4.1. Let D{x,y) be the formal power series defined by the equation 4>(:r, y; D(x, y)) = 0, 
where 


$(x,y; z) = z- 




/_ (y(l + xz) - (1 + y)z)(2 + xz) \ 
exp ^ x~ + xz ^ _(_ (i _|_ y)xz 2 )( 1 + xz) 2 ) ' 


Then, for every choice of y > 0 it holds that D(x,y) has a unique square-root singularity R(y) 
such that D{x , y) has a singular expansion of the following form in a domain dented at x = R(y): 


D(x, y ) = D 0 (y) + D 1 (y)X(y) + D 2 (y)X(y ) 2 + D 3 {y)X{yf + 0(X(y) 4 ), (13) 

where X(y) = y^l — x/R(y). In particular, for y = 1 we have the numerical values x = -R(l) = 
R « 0.05668, A,(l) « 1.82404 D^l) w -1.52769, D 2 (l) w 1.34779 and D 3 ( 1) « -1.25138. 


Proof. Fix y > 0. A simple computation shows that $ z (0, y; D(0, y)) = 1 > 0 and D(0,y) = y. 
Hence, by the Implicit Function Theorem, D{x,y) is analytic at x = 0. 

We start with showing that D(x,y) has a finite radius of convergence. Denote the singularity 
of the function D(x,y) by R(y). Observe that [x n ]D^(x,y) < [x n ]D{x,y), where Dq,(x,y) is the 
generating function associated with SP networks without a distinguished spanning tree. As it is 
shown in [S], the radius of convergence R%(y) of D${x,y) is finite. In particular, 0 < R(y) < 
R$(y) < 1 < oo and D(x,y) ceases to be analytic at x = R(y). 

Observe that the only source of singularity for D(x, y) is the condition 4 > z (R(y), y; D(R(y),y)) = 
0, meaning that the singularity arises from a branch point. Let us now justify that we have 
< f -zz(R(y),y, D{R{y),y)) 7 ^ 0 , which would give the claimed square-root expansion. This condition 
is enough in order to assure, for each choice of y , and square-root type singularity. For a contra¬ 
diction, let us assume the opposite. Hence, we have a solution (Rq, y 3 , zq) of the following system 
of equations: 

${x,y;z) = 0 , § z {x,y\z) = 0 , $ zz (x,y-,z) = 0 . 
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Observe that $(£, y; z) = z — A(x, y; z) exp(H(;r, y; z)) with A(x, y; z) and B[x , y; z) being rational 
functions. Hence, & z (x,y;z) = 1 — C(x, y; z) exp(B(a:, y; z)) where again C(x,y;z) is a rational 
function. Finally, & zz (x,y; z) can be written in the form E(x,y;z) exp(H(a;,y; z)) for a certain 
rational function E{x, y; z). 

In particular, combining the first two equations by eliminating the exponential term, we get 
the following system of rational equations: 

zC(x,y,z) = A(x,y;z), E(x,y,z) = 0. 

After rearranging the denominators in both expressions, such a system can be transformed into 
a system of two polynomial equations P\(x,y,z ) = 0, P 2 (x,y\z) = 0, from which we can get a 
new polynomial equation Q[x, y) = 0 by eliminating the variable 2 . By carrying out the explained 
computations with Maple, we obtain 


Q{x, y) = (-4y + yx- 4)y(y + 1 )T(x, y), 


where 


T(x, y) = 100(1 + y) 4 + 6917y(l + y) 3 x + 1266y 2 (l + y)V - 1867y 3 (l + y)x 3 + 280y V. 

We now argue that Q(x, y) = 0 does not have a solution with both y > 0 and x < 1. Observe that 
the first multiplicative term —4y + yx — 4 gives the solution (x, y) = (4 + 4/y, y). This means in 
particular that x is always greater than 1. It is also obvious that the multiplicative terms y and 
y + 1 cannot contribute with the required solution. Therefore, we need to analyse the existence of 
solutions (x,y) of T(x,y) with the condition y > 0 and a: < 1. Using that y(l + y ) 3 > y 3 (l+y), x > 
x 3 for all y > 0 and 0 < x < 1, we know that 6917y(l + y) 3 a: > 6917y 3 (l+ y)a ; 3 > 1867y 3 (l +y)x 3 . 
Hence, T(x, y) = 0 does not have solutions with both 0 < x < 1 and y > 0, which implies that the 
solution (xo,yo,zo) of the equation <&{x,y\z) = $ z (ir, y; z) = 0 satisfies that 4>, z (a;o, yo; -o) 7 ^ 0. 
Hence, the singularity of D{x,y) is of a square-root type in a domain dented at x = R(y). This 
proves the singular expansion in Equation 

In order to prove the special case of y = 1 in the statement of the lemma, we set y = 1, 
R = R( 1), and X = X(l) (and consequently x = R( 1 — X 2 )). By plugging the singular expansion 
of D(x , 1) in 4)(a:, y; z) = 0, taking the Taylor expansion in terms of X, and applying the method 
of indeterminate coefficients, we get the numerical values as claimed. Finally, observe that for 
each choice of y > 0, the generating function D{x,y) is aperiodic, as for every n there exists a 
network on n vertices. This gives that the singularity R{y) is unique, and the result holds. □ 


Knowing that D(x, y) admits a singular expansion of square root-type in a domain dented 
at x = R(y), one can compute by means of indeterminate coefficients the exact expressions of 
Di(y) for i > 1 in terms of the function D(R(y),y) = Z?o(y), which satisfies the functional 
equation $(i?(y),y, D 0 (y)) = 0. Although the expressions are long, we need to compute for 
enumerative purposes the evaluations at y = 1. For computational purposes, we include the 
following lemma where the coefficients of the singular expansions (rounded up to 5 digits) of 
the EGF D(x, 1), S(x, 1), 5(x, 1), P(x, 1), and P(x, 1) are obtained. Despite that in order to 
get asymptotic estimates for these counting formulas we only need the multiplicative constant of 
the term (1 — x/R) 1 ^ 2 , in order to get the precise asymptotic in the 2 -connected level we need 
expansions up to term (1 — x/R) 3 ^ 2 . 

Lemma 4.2. For each value of y > 0 the generating functions D, S , S, P, and P have a 
square-root singular expansion in a domain dented at R(y), where R{y) is the unique singularity 
of D(x,y). Furthermore, for y = 1 the singular expansions (with rounded coefficients) of D, S, S, 
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P, and P in a domain dented at R ss 0.05668 are 

D(x, 1) = A,(l) + D 0 (1)X + D 2 ( 1)X 2 + D 3 { 1)X 3 + 0(X 4 ) 
S(x, 1) = S 0 (l) + «Si(l)X + S 2 (1)X 2 + S 3 (1)X 3 + Q(X 4 ) 
S(x, 1) = S 0 (l) + £i(l)X + S 2 ( 1)X 2 + S 3 ( 1)X 3 + 0(X 4 ) 
P(x, 1) = P 0 ( 1) + Pi(l)X + P 2 (1)X 2 + P 3 (1)X 3 + 0(X 4 ) 
P(x, 1) = P 0 (l) + Pi(l)X + P 2 (1)X 2 + P 3 (1)X 3 + 0(X 4 ), 


where X = \/l — x/R, and the constants have the following approximate values: 



i = 0 

2 = 1 

i = 2 

i = 3 

A( 1) 

1.71871 

-1.17120 

1.17120 

-0.59820 

Si( i) 

0.17092 

-0.27289 

0.18433 

-0.15440 

Si( 1) 

0.30701 

-0.43079 

0.19616 

-0.12220 

m) 

0.65312 

-1.25480 

1.16347 

-1.09697 

Pi( i) 

0.41170 

-0.74041 

0.58941 

-0.47600 


Proof. The first claim follows directly due to Equations which are analytic and allow us 

to express D , S, S, P, and P in terms of D. In particular, all these generating functions have 
a unique singularity at x = R(y). The second part follows by setting y = 1 and by plugging the 
singular expansion of D(x , 1) into Equations (|6])-(11). □ 


Now we turn to the analysis of B(x,y), the EGF associated with the class of 2-connected SP 
graphs carrying a distinguished spanning tree. In our context, Equation § translates to 

2(1 + y)B v (x, y) = x 2 (l + D(x, y) + D(x, y) - y( exp (S(x, y) - 1))) (14) 

which roughly speaking means that when marking an edge in a 2-connected SP graph with a 
distinguished spanning tree, the resulting object is a network either of type D or P, but not a 
parallel network of type V with an edge linking the poles (see Equation (111). A direct integration 
of Equation (14) is technically involved due to the relations between the generating functions 
associated with the different types of networks. However, we can get a simple expression of 
B(x,y) in terms of the EGF associated with the networks just by combinatorial arguments. In 
the following lemma we provide such an equation. 


Lemma 4.3. The EGF B(x,y) associated with the class of 2-connected SP graphs carrying a 
distinguished spanning tree can be expressed as 

x 2 

B ( x , y ) = y v + Br ( x ' y ) + b m(x, y) - Br-m{ x : y), (15) 

where 

B R (x,y) = ^S(D-S), (16) 

2 

Bm(x, y) = {s( exp(S') - S - l) +yS(exp(S) - l) +y(exp(S) - S- l)), (17) 

B R _ M (x,y)= y(SP + SP). (18) 


Proof. Applying Tutte’s decomposition to 2-connected SP graphs bearing a distinguished spanning 
tree on at least 3 vertices only yields R-bricks (ring graphs) and M-bricks (multi-edge graphs), 
both carrying a distinguished spanning tree. In particular, there are no T-bricks since the set of 
h -networks is empty in our case. We obtain Expression (15) for B(x,y) using Equation © to 
which we needed to add x 2 y /2 since we also consider a single edge to be a 2-connected SP graph. 

Let us study each term. Let I? be a distinguished R-brick with a distinguished spanning tree. 
By definition, R is a cyclic chain of at least 3 networks that carries a spanning tree. In particular, 
exactly one of these networks is in V while the other ones are in V. This means that R can 
be decomposed into a non-series network in D and a series network in S that are joined by a 
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parallel operation and where the two poles are added to the graph, see also Figure [5] This gives 
Equation (16). 



Figure 5. Decomposition of a distinguished R-brick in the RMT-tree. 


We continue with M-bricks. Let M be a distinguished M-brick with a distinguished spanning 
tree. Then, M can be decomposed into at least three networks, all but possibly one of which are 
series and the possibly other one a single edge. These networks are joined by a parallel operation 
and the two poles are again added to the graph. This situation is similar to the decomposition 
of parallel networks carrying a spanning tree that we considered for developing Equation (10). 
The main difference is that, by definition, Ad is decomposed into at least three and not into at 
least two networks. We need to distinguish again between the two cases where there is a single 
edge component in M and where there is no such a component. Observe that it is not possible 
that there are two such components in M since we are only considering simple graphs. In the 
former of the two cases, we note that if the edge of the single edge component is contained in the 
distinguished spanning tree, then exactly one of the series networks is in S while all the others are 
in S. If, on the other hand, the edge is not in the spanning tree, then all series networks must be 
in S. This gives rise to Equation ( fl7| . 

Finally, we need to decompose 2-connected SP graphs with a distinguished spanning tree and 
with a distinguished {R, M}-edge in the RMT-tree. This means, that the distinguished edge 
corresponds to a virtual edge { x , y} matching a R-brick and a M-brick. Hence the graphs can be 
decomposed into a series network and a parallel network by a parallel operation, where we need to 
add again the two poles to the graph. One of the two networks must be in V while the other one 
must be in V. Hence, Equation (18) holds. See Figure [g] for an illustration of this situation. □ 




Figure 6. Decomposition of a distinguished {!?, M}-edge in the RMT-tree. 


We can now analyse the singular behaviour of B(x , y). 






SPANNING TREES IN RANDOM SERIES-PARALLEL GRAPHS 


13 


Lemma 4.4. Let y > 0. Then, B(x,y) has a unique square-root singularity, which is the unique 


singularity R(y) of the function D{x , y) in Lemma f.l Moreover, B{x , y) has a singular expansion 


of the following form in a domain dented at x = R{y): 


B(x,y ) = B 0 (y) + B 2 {y)X{yf + B 3 {y)X{y ) 3 + 0(X(y)% (19) 

where X(y) = \/l — x/R(y). When y = 1 we have the numerical values x = R{ 1) = R « 0.05668, 
B 0 (l) « 0.00176, B 2 ( 1) w -0.00394 and B 3 { 1) « 0.00062. 


Proof. Observe that the generating functions B^{x,y), B]\,.f(x,y) and B^_M{x,y) are analytic 
transformations of the generating functions for networks (namely, the EGFs in Lernma pUlj ). Hence, 
B(x,y) has a unique dominant singularity which is the same one as the coinciding singularity of 
EGFs of Lemma |4.2| namely R{y). Similarly, for each y, B(x,y) admits a singular expansion 
in a domain dented at R(y). In order to obtain it, we express the singular expansion of each of 
the network exponential generating functions appearing in Equation (15) in terms of the singular 
expansions obtained in Lemma |4.2| Observe that Equation (14) implies that the singular expansion 
of B(x, y) must start at X(y) 3 , which gives in particular that B\{y) = 0 (see Theorem VI.9 from 

ESI). 


Finally, by setting y = 1 and by the same procedure as above using Maple, we obtain the 
approximation of B,(\) for i > 0 as stated in the lemma. In particular, the term . 63 ( 1 ) depends 
on all singular coefficients in Lemma [472) □ 


Finally, we analyse the generating function C(x, y) of connected SP graphs carrying a dis¬ 
tinguished spanning tree. Since the singular expansion of B(x,y) is of a square-root type with 
exponent 1/2 as it is shown in Equation ( |19[ ), we get the singular expansion of C(x , y) immediately 
from Proposition 3.10 in [25] (see also |21j). 

Lemma 4.5. The singularity of C(x,y) is at p(y) = r(y)/exp(B x (T(y),y)), where r(y) is the 
unique solution of the equation T(y)B xx (r(y),y) = 1. The singular expansion of C(x,y) in a 
domain dented at p(y ) is 

C(x, y) = C 0 (y ) + C 2 {y)X{yf + C 3 (y)X(y ) 3 + 0(X(y)% 
where X ( y ) = \J\ — x/p(y) and 

Co{y) = t( 1 + \ogp(y) - log r(y)) + B{T(y),y), 

C 2 (y) = ~r(y), and 

c , x = 3 / 2p{y)exp(B x {p(y),y)) ~ 

3 V 2y TB xxx (r(y),y) - TB xx (r(y),y) 2 + 2B xx (r(y),y)' 

Additionally, when y = 1 we have p( 1) ~ 0.05288, Co(l) ~ 0.05450, C^l) = — r « —0.05668, 
andC 3 {l) « 0.00145. 

Proof. See the proof of Proposition 3.10 in ]25j for the analysis for a general value of y. When 
y = 1, we use Maple to obtain the approximations of the constants. The unicity of the singularity 
is assured by the aperiodicity of C{x , y) with y being fixed. See for instance the proof of [ 151 . 
Lemma 7] and [HI Lemma 9]. □ 


Now we have all necessary ingredients to prove the main theorem of this section. 


Proof of Theorem \1.1\ We will prove the statement for X n in detail. The result for Z n is obtained 
mutatis mutandis. Let us denote by C n the set of all connected SP graphs on n vertices, and C ^ 
the set of all connected SP graphs on n vertices carrying a distinguished spanning tree. For a 
graph G £ C n we write s(G) for the number of spanning trees in G. Then, the expected value of 
X n can be written as: 


E[X„] = Y, S ( G M G ] 

GeC„ 


^Gec„ S ( G ) 

\Cn\ 


\Cn\ 


[x n }C(x, 1) 

\Cn\ 


(20) 












SPANNING TREES IN RANDOM SERIES-PARALLEL GRAPHS 


14 


It follows directly from Lemma |4.5| and Theorem [37T] that the number of connected SP graphs on n 
vertices that carry a distinguished spanning tree is asymptotically equal to ^ 3 / 2 ^ n~ 5 ^ 2 p(l)~ n n\. 

The number of connected SP graphs on n vertices is asymptotically equal to c s n ~ 5 ^ 2 gj n n\, where 
c s ~ 0.0067912 and g s ~ 0.11021 are computable constants, as shown in [5( Theorem 3.7]. Di¬ 


viding the former by the latter as in Equation (20 1 , we obtain that the expected value of X n is 
asymptotically equal to sg~ n , where s « 0.09063 and g ^ 1 « 2.08415. 

The corresponding result for 2-connected SP graphs is obtained analogously by using Theo¬ 
rem 2.6 of g], which states that the number of 2-connected SP graphs on n vertices is asymptoti¬ 
cally equal to bn~ 5 ^ 2 r~ n n\, where b « 0.00101 and r ~ 0.12800. □ 

4.2. Fixing the edge density and limiting distributions. The previous results can be used to 
study random SP graphs with a fixed edge density, as well as limiting distributions for the number 
of edges. For sake of brevity, we only discuss the family of 2-connected SP graphs, but similar 
observations hold in the connected case. The first main important observation is that the number 
of edges in a uniformly at random 2-connected SP graph carrying a distinguished spanning tree 


follows a normal limiting distribution: Lemma 4.4 shows that the singular behaviour of B{x,y) 
is the same when choosing y in a real-valued neighbourhood of 1. Then, by the Quasi-Powers 
Theorem by Hwang (see |26j ) the distribution follows a normal limit law with linear expectation 
and variance. In particular, the number of edges is concentrated around its mean value. This 
behaviour is similar to the case of 2-connected SP graphs (without a spanning tree), where again 
the number of edges is normally distributed (see [5J). 

Under these circumstances, our techniques also provide a method to study the expected number 
of spanning trees in a random SP graph on n vertices of a given edge density g. Following the 
arguments of [Mi Theorem 3], for every g > 0 we can choose a value y 0 > 0 such that if we 
assign the weight i/q to each graph with k edges, then only the graphs with n vertices and with 
approximately gn edges (with a deviation of order n 1 / 2 ) have non-negligible weight. Such technique 
is valid whenever Quasi-Powers Theorem holds, and hence we can apply it in our context. 

As a case example, we plot the expected value of the random variable X U:fl that counts the 
number of spanning trees in a graph chosen uniformly at random from the 2-connected SP graphs 
with n vertices and edge density g. Let R(y) denote the radius of convergence of B(x,y). Given 
an edge density g, the right choice for yo is the unique positive solution of the following equation 
(see e.g. [MJ Theorem 3]): 

(21) 

Observe that when g tends to 1, the family of SP graphs under study are graphs with a small 
but positive number of cycles, whereas when g tends to 2 , the subfamily under study tends to the 
class of 2-trees. These cases correspond to the ones when y tends to 0 and infinity, respectively. 
Both cases will be analyzed fully in detail in Sections [5] and [6j respectively. 

The precise computational method to obtain the exponential growth constant of the expected 


value of X n ^ as a function of the edge density is the following. For a given density g we use (21 1 


to obtain the corresponding yo . Then we use the implicit expression of the singularity curve stated 
in Theorem 2.2. of 0 in order to obtain the growth constant of the number of 2-connected SP 
graphs of edge density equal to g. To get the growth constant in the setting of 2-connected SP 
graphs carrying a spanning tree and with edge density g , we perform the calculations explained 
in the proof of Lemma 4.1 for 3 / 0 - Finally, the exponential growth of the expected value of X ntfJj is 
obtained by dividing these two numerical values as we did in the proof of Theorem [Ll] In Figure [7] 
we plot this exponential growth constant in terms of the edge density g S (1.07626,1.97173). 

Let us mention that the non-plotted margins for g correspond to values of y very close to 0 and 
when y tends to infinity. In both cases, the numerical method used to get the constant growth 
for the number of spanning trees fails because of indetermination of the operation to be carried. 
Detailed analysis of the two cases when the edge density reaches its maximum and when it tends 
to its minimum will be carried out in Section [5] and in Section [6| respectively. 
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Figure 7. Exponential growth constant of the expected number of spanning trees 
in a random 2-connected SP graph (ordinate) as a function of its edge density 
(abscissa). 


4.3. Variance of the number of spanning trees in random SP graphs. Refining the com¬ 
binatorics exploited in the proof of Theorem one has also access to the second moment of 
the random variables X n and Z n . In this subsection we develop this by determining the growth 
constant of the variance of Z n . This will also show that Z n is not concentrated around its expected 
value. Recall that Z n was defined as the random variable which counts the number of spanning 
trees in a random 2-connected SP graph on n vertices. 

In order to determine the growth constant of the second moment of Z n , we will first study 
the asymptotic behaviour of the number of 2-connected SP graphs on n vertices carrying two 
distinguished spanning trees. As in Subsection |4.1[ we start with the analysis of networks carry¬ 
ing spanning trees and spanning forests. We define V* , S*, and V* as the classes of SP, series, 
and parallel networks, respectively, each carrying two distinguished spanning trees. Let D*(x,y ), 
S*(x,y), and P*(x,y) denote their EGFs. In order to be able to analyse these functions, we need 
again some auxiliary classes. Let V denote the class of all SP networks carrying a distinguished 
spanning tree and a distinguished spanning forest with (exactly) two components each of which 
contains one pole. Furthermore, let T> denote the class of all SP networks carrying two distin¬ 
guished spanning forests both with (exactly) two components each of which contains one pole. 
Let D(x,y) and D(x,y) denote their EGFs. In the same way S, V , S. and P as well as S(x,y), 
P(x,y), S(x,y ), and P(x,y ) are defined. We might again omit the parameters whenever they are 
clear from the context. Following the proof of Theorem 0 we start with the following lemma 
that provides the growth constant of these generating functions. 

Lemma 4.6. For y = 1, the generating functions D*, D, D, S*, S, S, P*, P, and P have a 
square-root expansion in a domain dented at R 2 « 0.02407. 

Proof. We start again with elaborating relations between all given generating functions. 

One can verify easily that the following relations hold: 


D*(x,y) 
D(x,y) 
D{x,y) 


y + S*(x, y) + P* (x, y), 
y + S(x, y) + P(x, y), 
y + S{x, y) + P(x, y). 


(22) 
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Recall that a series network can be decomposed into a network that is not series and an arbitrary 
SP network. Similarly to Equations ([8]) and we get the following equations for the generating 
functions associated with the networks from the classes S , S , and S. 


S*(x,y) = (D*(x,y) - S*(x,y))xD*(x,y), 

S{x, y) = (D(x, y) - S(x, y))xD*(x, y) + (D*(x, y ) - S*(x, y))xD(x, y), 

S{x,y) = (D*(x,y)~ S*(x,y))xD(x,y)+(D(x,y)~ S(x,y))xD*(x,y) 
+2 (D(x,y) - S(x,y))xD(x,y). 


(23) 


Finally, a parallel network can be decomposed into a set of at least one series network if the 
root edge is present, or into at least two series networks, otherwise. In the first case we need to 
distinguish whether the root edge is in a distinguished spanning tree or not. By a careful case 
distinction, we get the following three equations: 


P{x,y) = 


yS{x, y) exp(S(x, y)) 


P{x,y) 

P{x,y) 


S(x,y)(e-xp{S{x,y)) - l) + (1 + y)S(x, y) 2 exp(S(x, y)) 

+2yS(x, y) exp(S(x, y)) + y( exp(5(x, y) - l), 

{y + S{x, y)) ( exp(S(x, y)) - l) + yS(x, y) exp(5(a;, y)), (24) 

(exp(S(x,y)) - S(x,y) - l) +y(exp(S(x,y)) - l). 

Equations define a system of functional equations that can be analysed by means 

of Theorem 3.2 by setting y = 1. More precisely, each equation in this system is defined by an 
analytic function, because the exponential function is an entire function. In addition, easy lower 
and upper bounds imply that the radius of convergence of D(x, 1), D(x, 1), D(x, 1) is in (0, R], 
where R « 0.05668 is the constant obtained in Lemma |4~T) Solving the system of equations stated 
in Theorem |3.2| using Maple we get that f ?2 ~ 0.02407. Furthermore, Theorem 3.2 assures that 


i ?2 is the singularity of all generating functions appearing in Equations (22)-(24) and that they 


have a square-root expansion in a domain dented at R% ■ As all generating functions written so far 
are aperiodic, the singularity is unique on the circle |x| = i? 2 - □ 

Let B*(x,y) denote the EGF associated with the class of all 2-connected SP graphs carrying 


two spanning trees. As in Lemmas 4.3 and 4.4 one can write B*(x,y) as an analytic combination 
of the generating functions from Lemma 4.6 by using Tutte’s decomposition and Equation ([5]). 
Therefore, the dominant singularity of B* (x, 1) is the same as the coinciding singularity of the 
EGFs appearing in Lemma 4.6 which is R 2 , and moreover, B*(x, 1) has a square-root expansion 


in a domain dented at i? 2 - By Theorem |3.1| we get that the number of 2-connected SP graphs on 
n vertices carrying two spanning trees is asymptotically equal to 0 (n -5 / 2 I?^"”?r!). 

Let B n denote the number of 2-connected SP graphs on n vertices and SB* n the number of 
2-connected SP graphs on n vertices carrying two spanning trees. The second moment of Z n can 
be calculated in the following way: 

\SB* n \ [x n ]B*(x,l) 


E [Z 2 n ] = Y, s(G) 2 P[G] = 


GeB„ 


| B n 


| B n 


where again s(G) denotes the number of spanning trees in a graph G. Recall that the number 
of 2-connected SP graphs on n vertices is asymptotically equal to 0 (n -5 / 2 r -7l n!), where r « 
0.12800. This means that the second moment of Z n is asymptotically equal to 0 where 

w 2 1 ~ 5.31718. The same holds for Var (Z n ) = E [Z?] — E [Z n ] 2 since E [Z n ] is approximately 
equal to 0 (p _ra ), where p -1 ~ 2.08415 by Theorem |l.l| In particular, this implies that Z n is not 
concentrated around its expected value. 


5. Spanning trees in 2-trees 

In Subsection |4.2| we analysed the exponential growth constant of the expected number of 
spanning trees in a random 2-connected SP graph with respect to its edge density. Observe that if 
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the edge density of SP graphs tends to its maximum, we reach the class of 2-trees. In this section 
we present an alternative, direct way to compute this growth constant in the setting of 2-trees. 
More precisely, we give a proof of Theorem |1.2[ which is again based on the Symbolic Method, 
the extension of the Dissymmetry Theorem to tree-decomposable classes (Theorem |2.1[ ), and the 
singularity analysis of generating functions. As in the previous section, we use x and y to mark 
vertices and edges, respectively. In this particular scenario the number of edges is determined once 
fixing the number of vertices. However, for pedagogical reasons, we will write all the equations 
keeping track of both parameters. Recall that all graphs in this paper are considered to be labelled, 
unless otherwise specified. 

In order to obtain the expected number of spanning trees in a random 2-tree, we aim to 
determine the asymptotic precise estimate of the number of 2-trees as well as of the number of 2- 
trees carrying a distinguished spanning tree. As a first step, Lemma fO] provides the singularity qt 
and the singular expansion in a domain dented at qt of the generating function T(x : y) associated 
with the class of 2-trees. 

Lemma 5.1. For y = 1 it holds that T(x, y) has a unique square-root type singularity qt = 1/(2e) 
and admits the following singular expansion in a domain dented at x = qt : 

T(x, 1) = le- 3 / 2 - ^e- 3 / 2 X 2 + ^e" 3 / 2 * 3 + 0{X% 
where X = \J\ — x/ qt- 

Proof. Let T denote the class of labelled 2-trees rooted at an edge the endpoints of which do not 
bear a label and let T(x, y) denote its associated generating function. By the rules of the Symbolic 
Method for pointing operations we have the following relation between T(x,y) and T(x, y): 

y § y T ( x ’ y ') = y T ( x , y )- 

Observe that we had to add the factor x 2 /2 on the right-hand side because we had to add labels 
to the endpoints of the root edge. By integrating by substitution we get that 

T(*,„) = £ £ ^ (nx.y) - frfw) . 

In order to compute the radius of convergence of T(x,y ), it suffices to compute the radius of 
convergence of T{x,y) since their values coincide by the latter equation. A graph in T can be 
reconstructed by merging at the root edge a set of pairs of graphs in T that share a vertex, see 
Figure [S] Using the Symbolic Method this gives rise to the following equation for T(x,y): 

T(x, y) = j/exp (xT(x,y) 2 ). 

By Lagrange’s Theorem (see e.g. [53j) we get for every n,m> 0 that 



Figure 8. Decomposition of rooted 2-trees. 


m— li „mxu 


[x n }[y m ]T(x,y) = [x n ]-[u m - 1 ]e 


\u m 1 lm = x - 




0 


if n = 

otherwise. 
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We study the case y = 1 since the number of edges of a 2-tree is determined by its number 
of vertices. The inverse function of T(x, 1) is given by ip(u) = log [u)/u 2 . Let r > 0 such that 
?//(t) = (1 - 21og(r))/r 3 = 0, implying r = exp(l/2). Then we know from the Inverse Function 
Theorem that we have for the singularity qt of T(x, 1) that qt = V’( T ) = l/(2e). Therefore, 
T(x, 1) and hence also T(x, 1) have a square-root type singular expansion in a domain dented at 
qt- By the method of indeterminate coefficients we get that the singular expansion of T(x, 1) is of 
the form as stated in the lemma. Finally, by aperiodicity of the generating functions under study, 
the dominant singularity of T{x, 1) is unique. □ 


Now let us turn to 2-trees carrying a spanning tree. We denote this class by T s and its 
associated generating function by T s (x,y). Similarly to Section |4j we need to first analyse edge- 
maximal SP networks to get access to the singular behaviour of T s (x,y). For this purpose, let 
D r (x,y) denote the EGF associated with the set of all edge-maximal SP networks carrying a 
spanning tree that contains the root edge. Furthermore, let D r (x,y ) denote the EGF associated 
with the set of all edge-maximal SP networks carrying a spanning tree that does not contain the 
root edge. Finally, let D°{x,y) denote the EGF associated with the set of all edge-maximal SP 
networks with a distinguished spanning forest that consists of exactly two compon ents each of 
which contains one of the poles. Observe that we have D r (x,y) = D°(x,y). Lemma 
the singular behaviours of D r = D° and D r for y = 1. 


5.2 


gives us 


Lemma 5.2. We have D r = D° and for y = 1 the generating functions D r and D r have the same 
unique square-root singularity Rt ~ 0.07197. Furthermore, the singular expansions (with rounded 
coefficients) of D r [x, 1) and D r (x, 1) in a domain dented at x = Rt are: 

D r (x, 1) = 1.46516 - 0.77028X + 0.53282AT 2 - 0.40927X 3 + 0(X 4 ), 

D r (x, 1) = 0.34588 - 0.77028X + 0.8787CLY 2 - 0.92279X 3 + 0(X 4 ), 

where X = sj 1 — x/R t ■ 


Proof. Using the Symbolic Method (see also Figure [9]) one can easily verify that the following 
system of equations hold: 


D r (x,y) = yexp ( 2 (D r (x,y) + D r (x,y))xD r (x,y)^, 

D r {x, y) = yx(D r (x, y) + D r {x, y)) 2 exp (2 (D r (x, y) + D r (x , y))xD r {x , y)J . 




Figure 9. Possible decompositions of an edge-maximal SP network with a distin¬ 
guished spanning tree that contains the root edge (depicted on the left side)/does 
not contain the root edge (depicted on the right side). 


We may set y = 1 since the number of edges of 2-trees is always given by the number of 
vertices. The two equations in (25) are defined by entire functions. As these structures carry 
spanning trees, the corresponding singularity is smaller than qt, and hence their singularity Rt 


is finite. We can apply Theorem 3.2 and know therefore that D r (x, 1) and D r (x,l) have the 


same singularity Rt and that they have a square-root singular expansion in a domain dented at 
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this singularity. Solving the system of equations stated in this theorem with Maple yields that 
Rt ~ 0.07197. By aperiodicity of the counting formulas, this is the unique smallest dominant 
singularity. By means of indeterminate coefficients we get the exact coefficients of the singular 
expansions of D r (x , 1) and D r (x, 1) which are the ones as stated in the lemma. □ 

Lemma 5.3. Let y = 1. Then, T s (x, 1) has a square-root singularity, which is the singularity Rt 
of D r {x, 1) and D r (x , 1). Moreover, T s {x, 1) has a singular expansion of the following form in a 
domain dented at x = Rt ~ 0.07197: 

T s {x, 1) = Tq + T 2 s A' 2 + T 3 S X 3 + Q{ A' 4 ), 

where X = \J\ — x/Rt and Tq ss 0.00290, T| ss —0.00669, and T| ss 0.00133 are computable 
constants. 


Proof. We define the A e-tree t{G) of a graph G as the bipartite graph describing the incidences 
between edges and triangles of G. More precisely, the node set of r(G) is given by E(G ) U 
{{xi,x 2 ,x 3 } : {xi,Xj} G E{G) for each i 4 j G {1,2,3}} and two nodes e and A are neighbours 
if and only if e C A in G. 

Since we are dealing with graphs carrying a spanning tree we need to encode the information 
on the distinguished spanning tree also in the associated Ae-tree. For this reason, we define five 
different types of nodes of r(G) if G G T s ■ Let 14 and V e denote the set of vertices of r(G) 
associated with edges of G that are, respectively are not contained in the distinguished spanning 
tree of G. Moreover, for i G {1,2,3}, let Va, denote the sets of vertices of r(G) associated with 
triangles of G which contain exactly 3 — * edges that are in the distinguished spanning tree. 

Observe that for every G G T s the connected Ae-tree r(G) associated with it is uniquely 
defined. We claim that r(G) is a tree. Indeed, assume for a contradiction that there exists a 
cycle G = A 1 e 1 A 2 ... A k e k A 1 in r(G). Let G' be the induced subgraph of G on the vertex set 
|J ig r fc ] A 1 . In particular, being edge-maximal and A4-minor free, G' is also a 2-tree. However, G' 
does not contain a vertex of degree 2, a contradiction. As a consequence, r(G) is indeed a tree. 
This means that T s is a tree-decomposable class and therefore we can apply Theorem |2.1| Since 
Ae-trees are bipartite, the equation in Theorem 2.1 simplifies to T s — Tf — Tf_ 0 . The class Tf 

and 72 depending on whether the 


is naturally partitioned into the five classes Tf, Tf , T Al , T A2 
distinguished node is from the set 14, V e , V Al , Va 2 , or Va 3 . The class T 0 S _ 0 is partitioned into 
the classes Tf_ Ai , Tff_ A2 , 7}f_A, , 7{ s _a 2 > an d ^-a 3 by the structure of r(G). In terms of their 
associated generating functions these facts translate into 


T s =Tf + T f 


T: 


A, 


■T 


Ao 


rj-iS rjrS 

1 A 3 - 1 e-A 1 


- T£_ a 2 - t:_ Ai - T e _ a 2 - t:_ a 3 . (26) 

Using the Symbolic Method, it is not difficult to check that the following equations are true. Recall 
that we have D r (x,y) = D°(x,y). 


T£{x,y) = -B F {x,y), 


T s e =-D r {x,y), 


= y ( Dr ( x ’V)) 3 > 


T k 2 ( x >y) = x3 ( Dr i x ,y)) 2 Dr i x ,y), 


T As {x,y) = AD r ( x ,y ) (. D r {x,y)f 


and 


T i-A 1 ( x ,y) = x 3 (D r {x,y)) 3 , T£_ A2 (x,y ) = x 3 (D r {x, y)) 2 D r {x, y), 

T e-A 2 ( x ,y) = 2x3 (D r (x, y)) 2 D r (x,y), 


t La^V) = ^x 3 ( D r {x,y)) 3 


T e-A 3 (x,y ) = 2 x 3 D r {x,y)(D r (x,y)) . 

In particular, T s {x,y) can be expressed in terms of x, D r {x,y ), and D r (x,y) by plugging the 
previous equations in Equation (26). One can therefore easily verify that the dominant singularity 
of T s (x, 1) is the same as the coinciding one of D r (x, 1) and D r (x, y), namely Rt- Finally, we 
obtain the singular expansion of T s (x , 1) in a domain dented at R t by using Equation (26) and 
the singular expansions of D r (x, 1) and D r {x , 1) from Lemma 5.2 □ 


Now we have all ingredients that we need to finally prove Theorem 1 1.2 
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Proof of Theorem \l. 2\ Let U n denote the number of spanning trees in a random 2-tree on n 
vertices. Then, it holds that 


nu n ] = 51 s(G)P[G] 

Ger„ 


S 06 r, S (GQ 

\Tn\ 


I T°\ _ [x n ]T‘(x,l) 

\Tn\ \T n \ 


(27) 


where T n and denote the set of 2 -trees on n vertices and the set of 2 -trees on n vertices carrying a 
distinguished spanning tree, respectively, and s(G) denotes the number of spanning trees in a graph 
G. By Lemma [5.3| and Theorem |3.1| we get that the number of 2-trees on n vertices that carry a 
distinguished spanning tree is asymptotically equal to ^^ 3 / 9 ) n~ 5 ^ 2 R^ n n\. Furthermore, it follows 
from Lemma |5.f | and Theorem |3.f| that the number of 2-trees on n vertices is asymptotically equal 
to 4 ^ e _ 3 / 2 ) n~ 5 / 2 gy ra n!. Dividing the former by the latter as in Equation (27), we obtain that the 
expected value of U n is asymptotically equal to S 2 Q 2 n , where s 2 ~ 0.14307 and gf 1 ~ 2.55561. □ 


6. Spanning trees in series-parallel graphs with fixed excess 

Recall that the excess of a graph G is the number of its edges minus the number of its vertices. 
In this section we address the following question: given a fixed integer k (not depending on n), 
what is the expected number of spanning trees in a random connected SP graph with n vertices 
and excess equal to k if n is large? In order to study this question we exploit the structure of 
graphs with fixed excess introduced by Wright mm- The structure of graphs with fixed excess 
has been applied successfully in a wide variety of situations (see e.g. [3 m 131 [37]). 

We write Ck(x) for the EGF associated with connected SP graphs with excess equal to k and 
with a distinguished spanning tree. We denote by Ck{x ) the EGF of connected SP graphs with 
fixed excess equal to k. We also use the EGF for rooted labelled trees, that we denote by W(x). We 
recall that W( x) satisfies the functional equation W{x) = a; exp (IT (a:)). Concerning the singular 
expansion, W{x) has a unique square-root type singularity at x = e , and in a domain dented 
at this point W (x) has an expansion of the form 

W{ x) = 1 -V 2 Z + 0 (Z 2 ), 

where Z = \/l — ex. 

Let G be a graph with excess k. Starting from G, we build a multigraph that we call the kernel 
of G in the following way: we start deleting recursively vertices of degree one, thus obtaining the 
core of the graph (see for instance [34]). Then we continue by dissolving vertices of degree two, 
i.e. replacing the two edges incident to a vertex of degree two by a single edge. The resulting 
connected multigraph A'(G) (the kernel of G) has minimum degree greater or equal than 3 and 
fixed excess equal to k. The vertices of the kernel of G can be also labelled in the following way: 
the surviving \K(G)\ vertices in K(G) carry labels from 1 to n (if |G| = n ). The labels induce 
then an order of the vertices. The final labels in [1, |A'(G)|] of the kernel are the ones induced by 
these ordering of the primitive labels. 

If G has a distinguished spanning tree, the construction of the kernel of G induces also a 
spanning tree in K{G). Indeed, observe that the edges in G belonging to this distinguished 
structure induce a spanning tree in its core. Then, the distinguished spanning tree of the kernel 
is induced by the spanning tree in the corresponding core. In particular, possible loops do not 
belong to the induced spanning tree of the kernel. Additionally, if there is a multiple edge, at 
most one of the copies belongs to the spanning tree. See an example of the construction of the 
kernel in Figure flO] together with the distinguished spanning structure. 

It is straightforward how to reverse this construction: starting from a given kernel G we build 
directly the initial graph by pasting over each vertex of G a rooted labelled tree and substituting 
each edge by a sequence of rooted labelled trees. In our approach, dealing with spanning trees, we 
need to consider a slight modification for edges not belonging to the spanning tree. More precisely, 
we need to substitute in this case each edge (neither a loop nor a multiple edge) by 

1 + 2 W(x) + 3 W{xf + ■ ■ ■ = 5> + 1 )W{xY = (1 _^ (cc))2 , 
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Figure 10. Construction of the kernel of a graph. 


whose singular expansion in a domain dented at x = e 1 is of the form \Z 2 + 0(Z 4 ), where 
again Z = i/l — ex. The reason for this is that whenever we paste a sequence of rooted trees over 
the edge under consideration we must maintain the spanning structure acyclic and connected. If 
in the sequence we paste r different trees, the edge becomes subdivided into r + 1 edges. We need 
to choose then which one of these r + 1 edges is not in the spanning tree. 

In the case of the loop the previous sequence must start at the term 2W(x) in order to obtain 
at the end a simple graph. Similarly, for multiple edges we need to assure that at the end we 
finally obtain a simple graph. This case-by-case analysis suggests that generating functions would 
be somehow involved (see the similar problem in general graphs in mm)- However, we show 
that we are able to find closed formulas for the asymptotic estimates by means of a sequence of 
reductions. 

In the problem under study, our graph G is a connected SP graph, hence its kernel K(G) is 
a SP multigraph, i.e. a K 4 -minor free multigraph. It is also obvious that the kernel is a planar 
multigraph. Hence, by Euler’s relation, the number of possible kernels of excess k is finite. Observe 
that for a fixed value of k, a multigraph K of excess k maximizes its number of edges if and only 
if K is cubic (see the same argument in Lemma 4.3.1 of m for maps on surfaces). Hence, we can 
restrict ourselves to the study of connected SP graphs arising from a cubic kernel as these ones 
will provide the main contribution to the asymptotic. 

For technical reasons we will deal with weighted cubic multigraphs: the weight (called the 
compensation factor in |27], see also im has the following meaning: when substituting edges of 
the kernel by sequences of rooted trees, a loop has two possible orientations that give the same 
simple graph. A double (triple) edge can be permuted in two (six) ways and still produce the 
same object. The enumerative study of weighted cubic multigraphs will be postponed until the 
end of this section. We denote by gk the number of weighted cubic SP multigraphs with excess 
k. Additionally, we denote by g k the number of weighted cubic SP multigraphs with excess k 
that carry a distinguished spanning tree. In the rest of this section, all cubic multigraphs are 
considered to be weighted. 

Observe that a cubic multigraph of excess k has 2k vertices and 3 k edges. Therefore, any 
spanning tree has 2k — 1 edges and hence k + 1 edges are not used in the spanning structure. The 
first lemma gives a lower bound for \x n ]Ck{x): 


Lemma 6.1. For fixed k > 1, the following inequality holds: 

.(2,-2W{x)) k+1 W{x^ k 


g{k)[a 


(1 - WOe )) 4 *^ 1 


< [ x n ]C k {x ). 


(28) 


Proof. The proof is reminiscent to the proof of Lemma 3 in [35]. The EGF on the left hand side 
of Equation (28) can be written in the following way: 

, 2fe Wix )™- 1 


W( 3 


(1 - W(x)) 2k 


-1 


(3 — 2W (x)) k+1 W (x) 2k+2 
(1 - W(x)) 2k + 2 ' 


This can be interpreted as follows: for a given cubic SP multigraph with a distinguished spanning 
tree we 
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(a) paste a rooted labelled tree over each of the 2k vertices, 

(b) paste a sequence of at least one rooted labelled trees over each of the 2k— 1 edges belonging 
to the spanning tree of the kernel. 

(c) paste a sequence of at least one rooted labelled trees over each of the 2k + 2 edges not 
belonging to the spanning tree of the kernel, and then decide which of the new edges does 
not belong to the resulting spanning tree. 

Points (a) and (b) contribute with terms W(x) 2k and (W(x)/(1 — IE(:r))) 2fc_1 , respectively, where 
the second term is associated with a sequences of at least one rooted labelled tree. For the 
computation of Point (c), recall that loops are not in spanning trees. Hence, to guarantee that 
the final graph is simple we need sequences of length at least two for edges which do not belong 
to the spanning tree (length one is enough for multiple edges, but length two is needed for loops). 
Hence, the computation in point (c) arise from the fact that 

Y.<r + m(zY^ 3 - 2W(tm(x? 


r>2 


(1 -W{x)) 2 


This construction is injective and does not give all possible connected SP graphs with excess k. 
Summing all over all possible weighted cubic SP multigraphs we obtain the result as claimed. □ 

The next step now is to get an upper bound for \x n ]Ck{x). This is provided by the following 
lemma: 

Lemma 6.2. For fixed k > 1, the following inequality holds: 

W(x) 2k 


[.x n ]C k {x ) < g(k) [x 


(l + o(l)). 


(29) 


'{1-W(x)) 4k + 1 / 

Proof. The statement is proved by applying a similar argument to the one of Lemma |6.1[ For 
a fixed cubic multigraph with a distinguished spanning tree, we now paste over each edge an 
arbitrary sequence of rooted trees (possibly empty), and taking care of the special requirement on 
edges of the kernel such that do not belong to the initial spanning tree. In this way we generate all 
SP graphs carrying a spanning tree with excess k. Observe that it is possible to generate graphs 
which are not simple. Hence, this construction only gives an upper bound. Finally, the term o(l) 
in Equation (291 arises from the set of kernels that are not cubic. □ 


We can now combine both lemmas to get an asymptotic estimate for [x n ]Ck(x): 
Proposition 6.3. For fixed k > 1, the following asymptotic estimate in n holds: 

n 2k-l/ 2 


[x n ]C k (x ) = g k 


1 


2 2k+i/2 r(2fc +1/2) 


e"(l + o(l)). 


(30) 


Proof. We apply the Transfer Theorem for singularity analysis (Theorem 3.1l to Equations 
and (291. In particular we get that 


— [3 — 2W(x)) k+1 W(x) 6k 1 n 2k ~ 1/2 

*•'] Cdx)>^} { „ ( JL = 9 k 2 2fc + i/ 2 r(2 fc + 1 /2) 


(1 - W{x)) 4k+1 
and concerning the upper bound, 


HI)), 


W(a 


\2 k 


2k—1/2 


[x n ]C k (x) < g k {l* n ] {1 _ wix)) 4k + i ) C 1 + °(!)) = 9k 22k+1/ 2 r( 2 fc + 1 / 2 ) e ” (1 + ° (1)) ' 

In both estimates we have exploited the fact that W(e) = 1. As these estimates have the same 
singular behaviour we conclude the estimate in Equation (30). □ 


Before moving to the computation of ~g k an d let us mention that similar arguments give 
estimates for [x n ]Ck{x). Indeed, using the argument of Lemma 3 in [35] (or mutatis mutandis the 
previous arguments for C k (x)) one gets the following upper and lower bounds for \x n }Ck{x): 

W{xf k . r . A_m W[x) 2k 


9 k [X 


(1 - W{x)) 3k 


< [x n ]C k {x) < g k [x n ] 


(1 — W (x)) 3k 


(1 + 0 ( 1 )). 
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Again, by applying the Transfer Theorem for singularity analysis (Theorem |3.1| ) we get the estimate 

1 n 3k/2-l 

[x n ]C k (x) = g k 23fc/2 r(3fc/2) e"(l + o(l)). (31) 

Enumeration of cubic cores with spanning structures and proof of Theorem |1.3[ We 

complete the picture by getting asymptotic formulas (in k) for both g k and g k . In this section 
we get asymptotic estimates when k is large enough. The main results we implicitly use are 
the Transfer Theorems, joint with the fact that singular expansions can be integrated on dented 
domains (see for instance Theorem VI.9 in [23]). Additionally, given a value k , one can obtain 
the corresponding values by getting the Taylor expansions of the generating functions that will be 
introduced in the next lines. 

Let us start studying g k . As mentioned before in this section, g k is the number of connected 
cubic SP multigraphs with weights and excess equal to k. The weights are defined rigourously in 
the following way: given a multigraph with l-\ loops, I 2 double edges and I 3 triple edges, its weight 
is 2~ ll ~ l2 6~ l3 . Weights are needed to encode edge symmetries of multigraphs. Let G(u) be the 
EGF of connected weighted cubic SP multigraphs, where the variable u marks the excess. G(u) 
satisfies the following system of functional equations: 

= d(u) + c(u), 

'll 'll 

= -(d(u) + c(u)) + -, (32) 

= s(u) + p(u) + b(u), 
b(u ) 2 
u 

= c(u ) 2 — c(u)s(u), 

/ \ 1 / N 2 ^ 

= uc(u) + -uc(u ) + -. 

Full details concerning these equations can be found in Section 3 of [35], building on results 
on [28] (see also !)• The meaning of these EGF is the following: the term 6 u dG d ^ corresponds 
to the EGF of connected weighted cubic SP multigraphs where an edge (the root edge) is marked 
and oriented (remember that a cubic multigraph of excess k has 3 k edges and each edge has 2 
possible orientations). We have the following cases depending on the nature of the root edge: 
either the root edge is a loop (term b(u)) or a bridge (term d(u)) or when deleting it we get a 
series construction (term s(u)) or a parallel construction (term p(u)). The term c(u) plays the 
role of an auxiliary EGF. 

We proceed now to analyse this system. 


6 u 


dG(u) 
du 
b(u) 
c(u ) 
d(u) 
s(u) 
p{u) 


Proposition 6.4. The number g k = [u k ]G(u ) of weighted connected cubic SP multigraphs with 
fixed excess equal to k has the following asymptotic estimate: 


g k = ck 5/2 y fe (l + o(l)), 


where 7 = \J 6-\/3 — 9 ~ 0.17481 and c ~ 0.06034. 


Proof. From the system of equations (32) we get a single equation for c(it) by successive elimina¬ 
tion. Computations give that c(u) satisfies the algebraic equation 

0 = 8u + u 2 + (—8 + 24u + 6u 2 )c(u) + (—4 + 24u + 15u 2 )c(u) 2 + (8 u + 20 u 2 )c(u) 3 
+15 u 2 c(u) 4 + 6 u 2 c(u) 5 + u 2 c(u) 6 . 


We refer to Section VII .8 in [23]) for more details. From this equation we deduce that the dominant 
singularity of c(u) is the smallest positive root of the equation 

19683u 4 + 7776u 2 - 256 = 0, 

whose value is equal to 7 = ^\J 6\/3 — 9 « 0.17481 (—7 is also a root of this polynomial with 
the same absolute value, but it is easy to check that c(u) is analytic at u = —'})■ Using now 
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Newton’s Polygon Method (see Section VII.7 of [25] 1 we get that c{u) has a Puiseux’s expansion 
of the following form in a domain dented at u = 7 : 

c(u) = c 0 + CiU + 0(U 2 ), (33) 


where U = \J 1 — u/ 7 , cq « 0.61185 and Ci ~ —1.08766. Using Expansion (33) we deduce that 
D[u) admits the following singular expansion in a domain dented at u = 7 : 


d(u) = d 0 + diU + 0(U 2 ) 


with d 0 ~ 0.13306 and d\ ~ —0.19574. Finally, using that 6u d ° d ^ = d(u) + c{u) we conclude that 
the dominant singularity of G(u) is located at u = 7 . The proposition finally follows by integration 
of the Puiseux’s series (by applying Theorem VI.9 from .23]) in order to get the singular expansion 
of G(u) and by the application of the Transfer Theorem (Theorem |3.1[) □ 


We proceed with the study of g k . For this purpose, we refine the system of equations (32) 
in the following way: we denote by G(u) the EGF associated with connected weighted cubic SP 
multigraphs carrying a distinguished spanning tree (as before, u marks the excess). We study 

decompositions for 6 u dG ( j^f > , which correspond to the EGF of connected weighted cubic SP multi¬ 
graphs carrying a distinguished spanning tree where an edge is distinguished and oriented. 

In the following expressions, we use the subindex 0 to denote that the distinguished and oriented 
edge belongs to the spanning tree, while we use the subindex 1 to denote the opposite. Using this 
convention, we write do{u), bo(u), sq(u) and po(u) for the cases where this distinguished edge is 
a bridge, a loop, defines a series construction or a parallel construction, respectively, in such a 
way that the distinguished edge belongs to the spanning tree of the initial structure. Analogue 
definitions are done for di(rt), bi(u), si(u) and pi(u). Cq(u) and C\{u) are auxiliary families. 

Using the same arguments used to obtain the system of equations (32) yields the following, 
more involved system of equations: 


6u 


dG[u) 

du 

c 0 (u) 

ci(u) 

b 0 (u) 


bi{u) 


d 0 (u) 

di(u) 

s 0 (u) 

si{u) 

Po(u) 

Pi(u) 


d 0 (u) + c 0 (u) + di(u) + ci(u), 

b 0 (u) + s 0 (u) + po {u) , (34) 

bi(u) + si(u) +p 1 (u), 

0, 

U 'LL 

-(d 0 {u) + Co(u)) +u(di(u) + ci(rt)) + -, 
b\{u) 2 

1 

u 

0, 

(ci(u) - Si(u))ci(u) + ( c 0 (u ) - s 0 (u))c i(it) + (ci(u) - Si(u))c 0 (u), 

(ci{u) - Si(u))ci(u), 

'LL 1 

- + uc 0 (u) + 2uC\ (u) + -uc 0 (u) 2 , 
u + uco(u) + 3uci(u) + uco(u)ci(u) + 2uc\{u) 2 . 


We now analyse this system of equations similarly to what we did when studying System (32): 


Proposition 6.5. The number g k = [u k ]G(u) of weighted connected cubic SP multigraphs carrying 
a spanning tree and with excess equal to k has the following asymptotic estimate: 


g k = ck-^- k (l + o(l)), 


where 7 « 0.06709 and c « 0.06634. 
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Proof. The arguments are exactly the same as in Proposition |6.4| From the system of equa¬ 
tions (341 we get the following algebraic equation satisfied by co(u): 

0 = 121m 3 + 2304m + 7656m 2 

+(51696m - 26620m 2 - 968m 3 - 4608)c o (m) 

+(-384 + 34532m + 30888m 2 + 3388m 3 )c 0 (m) 2 
+(256 - 1392m - 10516m 2 - 6776m 3 )c 0 (m) 3 
+(32 - 1608m - 2376m 2 + 8470m 3 )c 0 (m) 4 
+(—352m - 132m 2 - 6776m 3 )c 0 (m) 5 
+(4m + 1144m 2 + 3388m 3 )c 0 (m) 6 
+(—44m 2 — 968m 3 )co(m) 7 + 121m 3 co(m) 8 . 

The singular point now is located at 7 w 0.06709, and the Puiseux’s expansion of Co(m) around 
m = 7 is of the form 

c 0 (m) = c 0 , 0 + c 0j i W + 0(W 2 ), 

where W = 1 — m/7, co,o ~ 0.29896, and co,i ~ —0.47032. We can then directly obtain the 
Puiseux’s expansion of Ci(m), &i(m) and di(u) from the expansion of co(m). By integra ting the 
equation 6ufj^G(u) = co(m) +ci(m) +do(w) and applying the Transfer Theorem (Theorem 3.1) we 
get the estimate for g k as it is claimed. □ 

We can now complete the proof of Theorem |1.3[ 

Proof of Theorem \1.3[ Due to Proposition 6.3 and Equation the value E[A'„,fc] is 


_ _ 2k —1/2 

[x n ]C k (x) _ 9k P( 2 fc+ i/ 2) era (1 + = g k P(3A:/2) 


[x n ]C k (x) 


rii 1 n 3k / 2 1 

yfc 2 3fc/2 r(3fe/2) e 


g k r(2fc + 1/2) 


©“( 1 + »(!))■ 


Using now the estimates obtained for g k and g k in Propositions | 6 .4] and 6.5 we obtain the result 
(and the constants) as claimed in the statement of the theorem. □ 


7. Concluding remarks and open problems 

In this paper, we have exploited the use of generating functions joint with analytic combinatorics 
to get exact expressions for the counting generating functions for connected SP graphs (2-connected 
SP graphs, 2-trees, and SP graphs with fixed excess) with a distinguished spanning tree. As a 
consequence, we were able to get very precise estimates for the expected number of spanning trees 
in an object chosen uniformly at random from the family under study. These techniques could 
be exploited in related families of graphs, as for instance the family of fc-trees with k > 3. Even 
though the analytic techniques used in this paper allow access to higher moments, unfortunately 
they do not provide a technique to determine the limit law of the number of spanning trees in SP 
graphs. 

With the techniques used in this paper one can also determine, with a bit more effort, the 
counting generating function of connected SP graphs carrying a distinguished spanning forest with 
a given number of components. From this one can derive the expected number of components of 
a random spanning forest in a random connected SP graph. In particular, the main difficulty in 
this situation is that for encoding networks carrying a distinguished spanning forest, one needs 
more auxiliary classes than in the restricted case of spanning trees. Roughly speaking, one needs 
to define classes of SP networks that carry a spanning forest and such that either the two poles 
are contained in the same component or in different components. The analysis of the generating 
function associated with networks as well as determining and analysing the generating functions 
associated with 2-connected and connected SP graphs carrying a spanning forest can be done 
similarly to the case of spanning trees. In the context of planar maps, let us mention that Bousquet- 
Melou and Courtiel m recently investigated the enumeration of regular planar maps carrying a 
distinguished spanning forest. 
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Finally, we would like to discuss similar results on planar graphs. Following the lines of [23] , 
the tools developed in the present paper can be extended easily to graphs defined by 3-connected 
components. When the family under consideration is defined by a finite number of 3-connected 
graphs, then the techniques used so far can be exploited to get analogue results. This would 
include, for instance, the family of graphs Ex(W 4 ) or Ex(W 5 ), where W 4 and W 5 are the wheel 
graphs with 4 and 5 external vertices, respectively (see [25]). In this context, a very interesting 
question is to extend the results of SP graphs to the random planar graph model. In this situation, 
the family of 3-connected components is infinite, and one requires extra results arising from map 
enumeration in order to control counting formulas for T-bricks. Although the number of maps 
carrying a spanning tree is a well-known fact (see e.g. El]), nothing is known when dealing with 
3-connected planar graphs. The problem of getting the expected number of spanning trees in a 
planar graph chosen uniformly at random will be considered in future investigations. 
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